Testing astrophysical models for the PAMELA positron excess with cosmic ray nuclei 
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The excess in the positron fraction reported by the PAMELA collaboration has been interpreted 
as due to annihilation or decay of dark matter in the Galaxy. More prosaically, it has been ascribed 
to direct production of positrons by nearby pulsars, or due to pion production during stochastic 
acceleration of hadronic cosmic rays in nearby sources. We point out that measurements of secondary 
nuclei produced by cosmic ray spallation can discriminate between these possibilities. New data on 
the titanium-to-iron ratio from the ATIC-2 experiment support the hadronic source model above and 
enable a prediction to be made for the boron-to-carbon ratio at energies above 100 GeV. Presently, 
all cosmic ray data are consistent with the positron excess being astrophysical in origin. 
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The PAMELA collaboration [l[ has reported an excess 
in the cosmic ray positron fraction, i.e. the ratio of the 
flux of positrons to the combined flux of positrons and 
electrons, 4> e +/{4>e+ + <fie-)i which is significantly above 
the background expected from production of positrons 
and electrons during propagation of cosmic ray protons 
and nuclei in the Galaxy [2j. It has been noted that 
the observed rise in the positron fraction between ~ 5 — 
100 GeV cannot be due to propagation effects [|[, rather 
it requires a local primary source of cosmic ray electrons 
and positrons, e.g. nearby pulsars 0, [E @, 0]- More 
excitingly, this could be the long sought for signature of 
the annihilation or deca y of dark matter particles in the 
Galaxy @, 1, M, EL dill ■ 

Alternatively the observed rise in the positron frac- 
tion could be due to the acceleration of positrons pro- 
duced by the decay of charged pions, which are created 
through hadronic interactions of cosmic r ay p rotons un- 
dergoing acceleration in a nearby source [14|. That the 
secondary-to-primary ratio should increase with energy 
if secondaries are accelerated in the same spatial region 
as the primaries has been noted quite some time ago in 
the context of cosmic ray acceleration in the interstellar 
medium [psl . [l6l . [ItJ ■ This model is conservative since it 
invokes only processes that are expected to occur in can- 
didate cosmic ray sources, in particular supernova rem- 
nants (SNR). One way to distinguish it from the other 
models is to e.g. compute the expected anti-proton-to- 
proton ratio, which is experimentally observed to be con- 
sistent with the standard background [18]. This is in fact 
in accord with the above model which predicts a rise in 
the p fraction only at energies above ~ 100 GeV [l9[ 
(see also ref.[2(|). This prediction cannot however be 
tested presently but must await data from the forthcom- 
ing AMS-2 mission [H| as well as PAMELA. 

Although dark matter annihilation or decay as the ex- 
planation of the positron signal would appear to be dis- 
favoured by the absence of a corresponding antiproton 
signal, this can in principle be accomodated in models 
with large dark matter particle masses or preferential lep- 
tonic annihilation/decay modes [12, [H, 0, [H| . Nearby 
pulsars as the source of the positrons are of course quite 
consistent with the absence of antiprotons. To differ- 



entiate between these possibilities and the model [14j in 
which secondary positrons from hadronic interactions are 
accelerated in the same region, we consider secondary nu- 
clei in cosmic rays which are produced by the spallation 
of the primaries. An increasing secondary-to-primary 
ratio (e.g. boron-to-carbon or titanium-to-iron) in the 
same energy region would confirm that there is indeed a 
nearby cosmic ray source where nuclei are being acceler- 
ated stochastically along wi th p rotons. 

An issue with this model [14| is that a crucial parame- 
ter is not known a priori but needs to be obtained from 
observations. This is the diffusion coefficient of relativis- 
tic particles near the accelerating SNR shock which de- 
termines the importance of a flatter spectral component 
over the usual Fermi spectrum and leads to the rise in 
secondary-to-primary ratios. Its absolute value cannot 
presently be reliably calculated. Observations of SNR in- 
dicate that the magnetic field is quite turbulent so that 
relativistic electrons diffuse close to the 'Bohm limit' with 
diffusion co-cfficicnt: £) Bohm = where the Larmor 

radius ri of the nucleus is proportional to the rigidity 
E/Z [1^|. We need to determine the actual diffusion 
co-efficient of ions in SNR in ratio to the Bohm value by 
fitting to data. A measurement of one nuclear secondary- 
to-primary ratio therefore allows us to make predictions 
for other ratios in the framework of this model. 

Very recently, data on the titanium-to-iron ratio 
(Ti/Fe) from the ATIC-2 experiment have been an- 
nounced [13| that indeed show a rise above ~ 100 GeV. 
We use this data as a calibration to determine the dif- 
fusion coefficient and then, extrapolating it according to 
its rigidity dependence, we predict the boron-to-carbon 
ratio (B/C) that should soon be measured by PAMELA 
(P. Picozza, private communication). 

Galactic cosmic rays with energies up to the 'knee' in 
the spectrum at ~ 3 x 10 15 eV are believed to be accel- 
erated by SNRs. The strong shocks present in these en- 
vironments allow for efficient diffusive shock acceleration 
(DSA) by the lst-order Fermi process (28[. In the sim- 
ple test-particle approximation, which is adequate for the 
level of accuracy of the present discussion [2{|, protons 
and nuclei that are injected upstream are accelerated to 
form a non-thermal power-law spectrum whose index de- 
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pends only on the parameters of the shock front, in par- 
ticular the compression ratio r. For a supersonic shock 
with r = 4 the steady state energy spectrum for protons 
and nuclei is NdE oc E t+2 dE where 7 = 3r/(r — 1). 
In the standard model of galactic cosmic ray origin, the 
accelerated primary nuclei produce secondaries by spal- 
lation on hydrogen and helium nuclei in the interstellar 
medium (ISM). In the simple 'leaky box model' [30( an 
energy-dependent escape of the cosmic rays out of the 
galaxy is invoked to obtain a secondary-to-primary ra- 
tio that decreases with energy as observed to date in the 
region ~ 1 — 100 GeV. This is also obtained by using 
the GALPROP code [2j which solves the full transport 
equation in 3 dimensions, and can yield both the time- 
independent as well as equilibrium solution. 

However, as the acceleration time for the highest en- 
ergy particles is of the same order as the timescale 
for spallation, the production of secondaries inside the 
sources must be taken into account. In any stochas- 
tic acceleration process one then expects the secondary- 
to-primary ratio to increase with energy since particles 
with higher energy have spent more time in the accelera- 
tion region and have therefore produced more secondaries 
[IM Il6l Il7j . This general argument can be quantified for 
the case of DSA by including the production of secon- 
daries due to spallation and decay as a source term, 
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S as ^ dec 



de fc , (1) 



where e k is the K.E./nucleon (in GeV), and a loss term, 
r ( N t (e k )de k = Nifa) K spall /3cn gas + — L_] de k , (2) 



where cr^f (cr| pa11 ) and rfe (r 4 dcc ) are the partial (total) 
cross-sections and decay time, respectively. The trans- 
port equation for any nuclear species i then reads 



Ox Ox o dx op 



(3) 



where fa is the phase space density and the different 
terms from left to right describe convection, spatial dif- 
fusion, adiabatic energy losses as well as losses and in- 
jection of particles from spallation or decay. We consider 
the acceleration of all species in the usual setup: in the 
frame of the shock front the plasma upstream (x < 0) 
and downstream (x > 0) is moving with velocity u_ and 
u+ respectively. We solve eq.Q analytically for relativis- 
tic energies e k greater than a few GeV/nucleon such that 
p rs E, /? » 1 and NidE ps 4np 2 fidp. At these energies 
ionization losses can be neglected and the spallation cross 
sections become energy independent. 

There are three relevant timescales in the problem: 



1 . Acceleration time r a 



U- — u + j Q \ u + u_ 



dp' 
~7 



(4) 



for Bohm diffusion and the parameter values men- 
tioned later. 

2. Spallation and decay time Tj. 



l/rf all ~1.2x 10 7 f-^) yr. (5) 



spall 



where an average cr, of 0(100) mb has been as- 
sumed. The rest lifetime r z doc of the isotopes con- 
sidered ranges between 4 x 10~ 2 yr and 10 17 yr. 

3. Age of the SNR under consideration [3, fl9| 

tsnr = x max /u + ~ 2 x 10 4 yr. (6) 

There are two essential requirements for SNR to effi- 
ciently accelerate nuclei by the DSA mechanism: 

(a) r acc <C T? pa , which is equivalent to 
20^^ < 1 =► £k «6.4x 10 5 ^ B^g GeV . (7) 



(b) tsnr *C n which implies, 
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« — => x — < 1 

I 11+ 



(8) 



The isotopes for which condition (b) is not satisfied at 
the lowest energy considered viz. 56 Ni, 57 Co, 55 Fe, 54 Mn, 
51 Cr, 49 V, 44 Ti and 7 Be do not contribute significantly, 
so their decays in the source region are neglected. 

We find that the general solution to eq. ([3]) for x ^ is 

ft = E (E%^ x/2 + F|e^/ 2 ) + Gf, (9) 
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where Xf = -^(l - y/l + 4Df rf/4), 
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where G i is the asymptotic value and Ejl and are 
determined by the recursive relations: 

a p± pit 

31 ~ D±\± 2 - 2u\f - 4r± ' 



F 



—4 v^ e^ 1 f^ 



(10) 

(11) 



We require that the phase space distribution function 
converges to the adopted primary composition Y, (at the 
injection energy po) far upstream of the SNR shock: 



fi(x,p) Yi6(p-p ) ,dfi/dp(x,p) x ^ °°> 0. 

We also require the solution to remain finite far down- 
stream. As the phase space density is continuous at the 
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shock front, we connect the solutions in both half planes 
to /P = fi(x = 0,p) and hnd them to be: 

3<i 

+ Y l 5(p-p )(l-e K < */ 2 ) , (13) 



ft = f?e xtx/2 + E E ii ( eAj+x/2 - cKx/2 



j<i 



+ G+(l-e A ^/ 2 ) 



(14) 



Using eqs.([7]|5]), we can linearly expand Xf and n i in 
eq. tfTO")) and the exponentials in eqs. (|13IU"4"|) 



3 A+x/2 ~ i _ ii_ e K ^/ 2 ~ (l + Li-V"-*/ 1 * (15) 



to obtain: 



/+ = /" + 



g+(a: = 0)-r+/« 



a; . 



(16) 



where g? 1 denotes the downstream/upstream source 
term: qf = V,. , fiH 

Finally we integrate the transport equation over an 
infinitesimal interval around the shock, assuming that 
(it lii = r t/ r 7 = «ga S /™gas = r and that D+ ~ 

^ = -7/f-7(l + r- 2 )^/f 



+ 7 



(l + r 2 )^fcffi 



•5^^(p-po) 



(17) 



which is readily solved by 



/0(p) = I ?LY e-T( ] -'' ,U) . : " 1,1 



o P \ P 



x 7 



gr(x = 0)£) i -(p') 



+ y^(p'- Po ) 



(18) 



Our eqs. (|16H18"|) should be compared to eqs. (4-6) of 
ref.[l4[ where the loss terms Tifi were not taken into ac- 
count. The exponential in our eq. (fTB)) leads to a natural 
cut-off in both the primary and secondary spectra above 
the energy predicted by eq. ([7]). However, due to the 
approximations we have made, the secondary-to-primary 
ratios cannot be predicted reliably for ATiDi/u 2 > 0.1 
i.e. much beyond ~ 1 TeV. 

Starting from the heaviest isotope, eqs. lfll))) and (fT5|) 
can be solved iteratively to obtain the injection spectrum 
after integrating over the SNR volume, 

Ni{E)=4n dxp 2 f i (p)4 7 rx 2 . (19) 

Jo 

To account for the subsequent propagation of the nuclei 
through the ISM we solve the transport equation in the 



'leaky box model' [30j which reproduces the observed de- 
crease of secondary-to-primary ratios with energy in the 
range ~ 1 — 100 GeV by assuming an energy-dependent 
lifetime for escape from the Galaxy. The steady state 
cosmic ray densities A/i observed at Earth are then given 
by recursion, starting from the heaviest isotope, 



Mi = 



l/r e! 



(20) 



where TZsn ~ 0.03 yr -1 is the Galactic supernova rate. 

We calculate the source densities iVj and ambient den- 
sities Afi, taking into account all stable and metastable 
isotopes from 64 Ni down to 46 Cr/ 46 Ca for the titanium- 
to-iron ratio, and from 18 down to 10 Be for the boron- 
to-carbon ratio. Short lived isotopes that ^ decay im- 
mediately into (meta)stable elements are accounted for 
in the cross-sections. The primary source abundances 
are taken from ref.[3l| and we have adopted an injection 
energy of 1 GeV independent of the species. The partial 
spallation cross-sections are from semi-analytical tabula- 
tions [32| and the total inelastic cross-sections is obtained 
from an empirical formula [33j . The escape time is mod- 
elled according to the usual relation: 



pcx c 



PCXescA E / Z i 



(21) 



where x is the column density traversed in the ISM and 
p = 0.02 atom cm -3 is the typical mass density of hydro- 
gen in the ISM. We have neglected spallation on helium 
at this level of precision as its inclusion will have an ef- 
fect < 10%. The parameters of the fit are sensitive to 
the adopted partial spallation cross-sections, for exam- 
ple n ~ 0.7 for the Ti/Fe ratio but - 0.6 for the B/C 
ratio. For improved accuracy a current compilation of 
experimentally deduced cross-sections should be used. 
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FIG. 1: The titanium-to-iron ratio in cosmic rays along with 
model predictions — the 'leaky box' model with production 
of secondaries during propagation only (dashed line) , and in- 
cluding production and acceleration of secondaries in a nearby 
source (solid line - dotted beyond the validity of our calcula- 
tion). The data points are from ATIC-2 (triangles) [27t ] and 
HEAO-3-C3 (circles) @. 
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FIG. 2: The boron-to-carbon ratio in cosmic rays along with 
model predictions — the 'leaky box' model with production 
of secondaries during propagation only (dashed line), and in- 
cluding production and acceleration of secondaries in a nearby 
source (solid line). The data points are from HEAO-3-C2 (cir- 
cles) HI, ATIC-2 (triangles) [H and CREAM (squares) [H. 

Following ref . (l9| , the parameters are chosen to be: r = 
4, u- = 0.5 x 10 8 cms _1 , n^ s = 2cm~ 3 and B = 1 /iG. 
The diffusion coefficient in the SNR is 

Di(E) = 3.3 x lO 22 ^- 1 B^EdevZr 1 cmV 1 (22) 

where the fudge factor is the ratio of the diffusion 
coefficient to the Bohm value and is determined by fitting 
to the measured titanium-to-iron ratio. 

The calculated titanium-to-iron ratio together with rel- 
evant experimental data is shown in Fig. [T] The dashed 



line corresponds to the leaky box model with production 
of secondaries during propagation only and is a good fit to 
the (reanalysed) HEAO-3-C3 data [34[. The solid line in- 
cludes production and acceleration of secondaries inside 
the source regions which results in an increasing ratio 
for energies above ~ 50GeV/?i and reproduces well the 
ATIC-2 data [27| taking^" 1 ~ 40. This is similar to the 
value reported in refs.jlil Il9j. thus ensuring consistency 
with the e + as well as p fraction measured by PAMELA. 

Clearly the experimental situation is inconclusive so 
a new test is called for. Fig. [3] shows the correspond- 
ing expectation for the boron-to-carbon ratio with the 
diffusion coefficient scaled proportional to rigidity ac- 
cording to ea.(|22"l). The CREAM data [3gl do show a 
downward trend as has been emphasized recently [37l ]. 
but the uncertainties are still large so we await more 
precise measurements by PAMELA which has been di- 
rectly calibrated in a test beam [38|. Agreement with 
our prediction would confirm the astrophysical origin of 
the positron excess as proposed in ref.[14| and thus es- 
tablish the existence of an accelerator of hadronic cosmic 
rays within a few kpc. 
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